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We present results for QCD with Nf — 2 flavors of dynamical quarks using nonperturbatively improved Wilson 
fermions at finite temperature on 16 3 x 8 and 24 3 x 10 lattices. We determine the transition temperature in 
the range of quark masses 0.6 < m w /m p < 0.8. After fixing the Maximal Abelian gauge we investigate the 
contribution of Abelian monopoles to the Polyakov loop, Polyakov loop susceptibility and confirm Abelian and 
monopole dominance in full QCD. 



1. INTRODUCTION 

In order to obtain predictions for the real world 
from lattice QCD, we have to extrapolate the lat- 
tice data to the continuum and to the chiral lim- 
its. Recently the Bielefeld group and the CP- 
PACS collaboration [2] using different fermion 
actions obtained consistent values for the criti- 
cal temperature T c in the chiral limit, albeit on 
rather coarse lattices at N t = 4 and 6. Edwards 
and Heller [3] determined T c for N t = 4, 6 us- 
ing nonperturbatively improved Wilson fermions. 
We compute T c on finer lattices with N t = 8 and 
10 with high statistics. Our results for N t = 8 
were reported in Ref. 0]. 

2. SIMULATION 

We use fermionic action for nonperturbatively 
0{a) improved Wilson fermions: 

Sf = Sp- ^Kgc sw a 5 ^2^'ip{x)cr l _ lu F llu ip(x), 

where Sp is the original Wilson action, c sw was 
calculated in 0. 



Configurations are generated on 16 3 x 8 (f3 — 
5.2 and 5.25) and 24 3 x 10 (/3 = 5.2) lattices at 
various k. The values of n and corresponding 
number of trajectories for 16 3 x 8 and 24 3 x 10 lat- 
tices can be found in Ref. 0] and Table 1, respec- 
tively. The number of configurations for 24 3 x 10 
lattice is not large enough and results for this lat- 
tice are preliminary. 



K 


0.1352 


0.1354 


0.1356 


0.1358 


# traj. 


5,400 


7,400 


3, 130 


1,650 



Table l:Simulation statistics on 24 3 x 10. 
We use results obtained at T=0 to fix the scale. 
An updated version of the contour plot of lines of 
constant r^/a and m^jm p |5] is shown in Fig.^ 

3. CRITICAL TEMPERATURE 

We use non-Abelian and Abelian Polyakov 
loops to determine the transition temperature. 
The Maximally Abelian gauge [7] is fixed by max- 
imizing the quantity 8 
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"Talk given by Y. N. at Lattice'03. 



with respect to gauge transformations g, 
U(s,fi) — g(s)U(s, /i)5 t (s -I- jl). We use an 
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Figure 1. The lines of constant r$/a and m^/mp 
at T — 0. Crosses correspond to parameters used 
in this work. 

SU(3) version of the simulated annealing algo- 
rithm [5]. The Abelian link variables are defined 
as Ui(s,iJ,) = exp{i#* (s)}, where 

6^(s) = arg Un(s,p) - -)] arg(7 j j(s,^)| mo d2,r • 

We then define the Abelian, monopole and pho- 
ton Polyakov loop operators as in ^U] . 

The Polyakov loop susceptibility is used to de- 
termine the transition point. In Fig[3the non- 
Abelian Polyakov loop and its susceptibility are 
depicted. We get the following values for the crit- 
ical temperature: 



T c ~ 196(8)MeV, 
T c = 210(4)MeV, 
T c = 219(3)MeV, 

It is known JI] 



m^/mp r 
m^/nip 



0.64 (Preliminary) 
= 0.77 
= 0.81 



that in SU(2) gluodynamics 
the non- Abelian, Abelian and monopole Polyakov 
loops give the same T c and critical indices. It fol- 
lows from comparison of Fig[3and Fig[3|that also 
in the full QCD these three Polyakov loops have 
similar behavior and their susceptibilities have 
maxima at the same temperature. 

4. CONTINUUM LIMIT 

At small enough lattice spacing and quark mass 
one can extrapolate the critical temperature T c to 
the continuum and the chiral limits using formula: 

T c r = (T c r ) m ^° + C a (a/r ) 2 + C q (±- - 
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Figure 2. Non- Abelian Polyakov loop (top) and 
its susceptibility (bottom). 
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Figure 3. The same as in Fig. but for Abelian, 
monopole and photon contributions. 

where ro = 0.5 fm and (T c ro) m<!,a ^ corresponds 
to the extrapolated value. 

We are brave enough to use four values for 
T c ro (see Table 3), obtained at rather large quark 
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masses, to estimate the parameters in this ex- 
trapolation expression. A fit gives (T c ) m <" a ^ 
and a with large errors: (T c ) m i' a ^° ~ 190 MeV, 
a - 0.8. 



T c r 


a/r 




0.50(2) 


0.20(1) 


N t = 10, (3 = 5.2 (prelim.) 


0.53(1) 


0.234(5) 


N t =8, (3 = 5.2 


0.55(1) 


0.225(5) 


N t =8, f3 = 5.25 


0.57(2) 


0.29(1) 


N t = 6, = 5.2 (Ref. EJ) 



Table 3: Available data for T c rQ. 

In Nf — 2 QCD the critical indices are ex- 
pected to belong to the universality class of the 
3D 0(4) spin model for which one expects a = 
0.55. Fixing a to this value we get the extrapo- 
lated temperature with a higher accuracy: 



r ™„a^o _ 172.5(3.3) MeV. 



(1) 



This value agrees with values obtained in Refs. ^ 

m 

5. CONCLUSIONS 

We determined the critical temperature in full 
QCD on 16 3 x 8 and 24 3 x 10 lattices at /3 = 5.2 
and 5.25 with Nf = 2 clover fermions using non- 
Abelian and Abelian Polyakov loop susceptibili- 
ties. Our results are in agreement with the results 
of other groups, as it is shown in Fig. Assum- 
ing that the critical indices of the two flavor QCD 
belong to the universality class of the 3D 0(4) 
spin model, we extrapolate the critical tempera- 
ture T c to the continuum and the chiral limits. 
The extrapolation result is given by Q. We are 
continuing simulations on 24 3 x 10 lattice to get 
better precision of T c value on this lattice. 
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Figure 4. T c J^fa, as a function of m^/rrip. The 
circle and squares show results of an d re ~ 
spectively. The diamond and triangles correspond 
to our data. 
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